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Abstract 

We have analyzed the joint charged particle to photon density correla- 
tions in WA98 experiment at CERN-SPS in order to study the formation 
of the chirally disoriented vacuum in Pb + Pb collisions at 158 AGeV. 
The method of analysis is based on the event-by-event estimation of a set 
of robust observables which are the ratios of normalized bivariate facto- 
rial moments corrected for many efficiency factors. This is found to be 
a very powerful technique in identifying, within the limits of the detector 
resolution, the events with the localized charged particle to photon density 
fluctuations. 



In high energy heavy ion collisions, hot fireballs may be produced with in- 
teriors having anomalous chiral order parameters. Eventually, these conjectured 
domains of disoriented chiral condensate(DCC) would decay by radiating pions 
with distinctive isospin distribution, following an inverse square-root law |l| 
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where / is the neutral pion fraction, given by 
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This is markedly different from the generic distribution, which is governed by a 
binomial distribution sharply peaked at / ~ 1/3. 

Therefore, an experimental study of event by event charged to neutral density 
correlation would shed light on this exotic phenomenon. 

WA980 is the only heavy ion experiment at CERN SPS which has initiated 
the 'vacuum engineering' in Pb + Pb collisions at 158 AGeV. The multiplicities of 
the charged particles are measured by a Silicon Pad Multiplicity Detector(SPMD) 
within a pseudorapidity coverage of 2.35 < rj < 3.75. Photon multiplicities are 
measured in the range of 2.4 < 77 < 4.4 by a highly granular Photon Multiplicity 
Detector (PMD). The efficiency of detecting a charged particle in the active area of 
the detector is ~ 99%, while the average photon counting efficiency is 70% — 80% 
depending on centrality, with a 35% contamination of showering hadrons. The 
centrality of the collision is determined by a Mid Rapidity Calorimeter (MIRAC) 
and a Zero Degree Calorimeter(ZDC) by respective measurement of the transverse 
and forward energies. 

In addition to the detection efficiencies for charged particles and photon mea- 
surements, the following probability factors also exist because of finite acceptance 
of the detectors: 

(a) e c h- Probability of observing a n ± within the detector acceptance, (b) 
eo: Probability that none of the photons from a 7r° falls in the detector, (c) t\\ 
Probability that one of the 7's from a it comes within the detector coverage, (d) 
t2- Probability of observing both the 7's from a 7r° within the detector acceptance. 

These e factors along with the detection efficiencies and the background 
effects0 put some serious challenges for us in this game. 

Nevertheless, a sophisticated analysis technique, based on the method of the 
normalized factorial moments, has been developed to estimate, within limits, 
the inclusive charged to neutral density correlations. In this formalism, the ac- 
ceptance factors, e c h and (i=0,l,2), which can not be directly estimated from 
the data, are incorporated in the generating function of the bi-variate factorial 
moments. The resulting factorial moments contain the efficiency factors. 

It can be shown 0] that the proper combinations of these factorial moments 
eliminate the acceptance factors and formulate the following 'robust' observables 



In the equations cited above, the angular brackets < ... > denote the averages 
over a large number of events and N is the individual event multiplicity. 
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where F, L and F^j are given by: 
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Moreover, one finds that for alH > 1 
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in the two cases. 

The above formalism developed by the MiniMax Collaboration is not sensitive 
to the localized fluctuations in single events. In order to make it suitable to 
identify the single event fluctuation structures, we have to calculate the robust 
observables on an event-by-event basis. In that case the angular brackets in the 
numerators and denominators of equations (4) - (5) should denote the averages 
over a considerable number of bins in a single event phase space distributions 
of N 1 — N c h. Now, the question comes whether one can use equation (6) to 
distinguish between a pure DCC and pure Generic event. The issue here is the 
size and correlation length of a DCC domain |J. For example, if there is a single 
domain covering the entire detector acceptance, then one would expect that each 
bin (up to statistical fluctuations) would have the same orientation. On the other 
hand, if the DCC were produced throughout the acceptance but the correlation 
length is comparable to the bin size, then one might expect a l/y/f distribution 
(contaminated by whatever background is present) more or less bin by bin and 
in that case equation (6) can be used to distinguish a DCC and a generic event. 
The purpose of the multi-resolution analysis with moments or wavelets is to find 
out this type of correlations. Hence, in the WA98 environment we shall compute 
the robust observables event - by - event and try to explore their properties. This 
is important for using them as a distinctive tool for a DCC and generic event. 

To analyze the data, we have selected the central events which correspond to 
about top 10% of the minimum bias cross-section. The selection of the central 
events confine us to the high detection efficiency regions specially for photons. 
The joint charged particle / photon distributions were constructed by using the 
overlap region 3.25 < rj < 3.75 of PMD and SPMD for the full azimuthal coverage. 
We have partitioned the whole azimuthal coverage into four different types of bins, 
namely, (a) 32 bins: A0 = 11.25°, (b) 16 bins: A0 = 22.50°, (c) 8 bins: A0 = 
45.00°, (d) 4 bins: A0 = 90.00°, in order to achieve the multi-resolution binning 
of the available i] — <p region. 

In order to search for any anomalous charged to neutral particle density fluc- 
tuation, we have computed the robust observables for the following three sets of 
events: (1) generic events, (2) simulated exotic events and (3) experimental data. 

The generic events are generated by the Monte-Carlo code VENUS 0. The 
detector effects are introduced in VENUS in order to generate data-like generic 
events. 

The simulation of the exotic events has been done by modifying the output of 
the VENUS event generator to include the fluctuations in the relative production 
of the charged particles and neutral pions. We assume that a single DCC domain 
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translates to the percent of VENUS pions that gets affected. The identities of 
these affected pions are altered by choosing a certain value of / randomly from 
the distribution given by (f). For the present purpose, we have simulated the 
exotic events in such a way that all of the pions within 3.0 < rj < 4.0 and 
0° < <fi < 45° are supposed to come from a DCC domain. This is an example 
of a very large domain of DCC which may not be created in experiments but is 
very useful to test the sensitivity of the method towards the localized charged to 
neutral fluctuations. 

The results of the analysis are summerized in Table -1 for the above three 
cases. The entire analysis has been done with a large number of events so that 
the statistical errors are small. In case of the data and simulated generic events, 
the distributions of r^i are peaked around 1.00. But the mean value of the 
distributions shift to a lower value in case of the simulated exotic events. This 
indicates the sensitivity of the method to disentangle exotic events accompanied 
with a large domain of DCC from the generic samples. Moreover, the difference of 
the data and the simulated exotic events indicates the probability of the formation 
of smaller and/or multiple domains of DCC. Further simulations in this direction 
are in progress. The sensitiveness of the higher order moments (say, r 31 ) to the 
localized fluctuations is remarkable within the limits of the so called 'empty bin 
effects'. The root mean square values of the distribution of r^i, for all resolutions, 
are higher for data than the simulated normal events. Also, the values of a 
increases with decreasing bin size. This indicates that one should not go up to 
the smallest bin samples to avoid the trivial statistical effects. 

Table - 1 

The average values of r iti for the simulated normal events, simulated exotic 
events and data. 



Event type 


INo. of Bins 


H,l ± cr 


r 2 ,i ± a 


7*3,1 ± cr 


JNormai 

Exotic 

Data 


4 


1.001 ±0.02 
0.845 ± 0.04 
1.009 ± 0.02 


0.997 ±0.02 
0.806 ± 0.04 
1.019 ±0.03 


1.008 ±0.03 
0.742 ± 0.05 
1.031 ±0.05 


JNormai 

Exotic 

Data 


8 


1.001 ±0.01 
0.852 ± 0.04 
1.012 ±0.03 


1.004 ±0.02 
0.848 ± 0.04 
1.027 ±0.05 


1.009 ±0.04 
0.801 ±0.05 
1.044 ±0.07 


JNormai 

Exotic 

Data 


16 


1.002 ± 0.02 
0.889 ± 0.04 
1.021 ±0.05 


1.007 ±0.04 
0.853 ± 0.04 
1.046 ± 0.07 


1.010 ±0.06 
0.831 ±0.05 
1.080 ±0.11 


JNormai 

Exotic 

Data 


32 


1.002 ± 0.04 
0.901 ± 0.04 
1.028 ±0.05 


1.009 ± 0.08 
0.889 ± 0.05 
1.067 ±0.11 


1.090 ±0.12 
0.821 ±0.08 
1.130 ±0.19 



A closer look into the data indicates a few events having the values of r i)X > 3<x, 
where a denotes the root mean square value of the concerned distribution. These 
events may be interesting with the indication of considerable charged particle 
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/ photon density fluctuations in the azimuthal space. But a concrete inference 
needs more studies on the detector effects as the data of both the PMD and 
SPMD are in a preliminary stage. This is, therefore, evident that in the first 
attempt of event-by-event computation of the robust observables, we have estab- 
lished the suitability of the moment analysis in the DCC search. It is found that 
the data contains more fluctuations than the simulation and for both the data 
and simulation a increases with decreasing bin size. Also, the preliminary data 
indicates the presence of remarkable charged to neutral density fluctuations in 
some events. These results, thus, demand a more careful analysis of the detector 
effects for a confident conclusion. 
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